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ONE-DIMENSIONAL MODEL EQUATIONS FOR HYPERBOLIC 

FLUID FLOW 

TAM DO, VU HOANG, MARIA RADOSZ, AND XIAOQIAN XU 


Abstract. In this paper we study the singularity formation for two nonlocal ID active 
scalar equations, focusing on the hyperbolic flow scenario. Those ID equations can be 
regarded as simplified models of some 2D fluid equations. 


1. Introduction 
The following transport equation 

(1) uj t + u ■ Vca = 0. 

is a basic mathematical model in fluid dynamics. If u depends onw, fll]) is called an active 
scalar equation. The problem of deciding whether blowup can occur for smooth initial 
data becomes very hard if the dependence of u is nonlocal in space. 

The relationship expressing u in terms of uj is commonly called Biot-Savart law. We 
have the following examples in 2D: 

(2) u = V ± (-A)" 1 cu, 

where V x = (— d y: d x ) is the perpendicular gradient. Equations (JTj) and (J2]) are the 
vorticity form of 2D Euler equation. When we take 

u = V ± (— 

(JTJ) becomes the surface quasi-geostrophic (SQG) equation, which has important applica¬ 
tions in geophysics, or can be regarded as a toy model for the 3D-Euler equations. For 
more details we refer to [I]. 

A question of great importance is whether solutions for these equations form singular¬ 
ities in finite time. A promising new approach for the construction of singular solutions 
is to use the hyperbolic flow scenario. In 0, 0. such a scenario was proposed to obtain 
singular solutions for the 3D Euler equations, and in [10], the long-standing question of 
existence of solutions to the 2D Euler with double-exponential gradient growth was settled 
using hyperbolic flow. 

The hyperbolic flow scenario in two dimensions can be explained in the following way. 
Consider e.g. a flow in the upper half-plane {X 2 > 0}. The essential properties required 
are (see Figured] for an illustration): 

• There is a stagnant point of the flow at one boundary point (e.g. the origin) for 
all times. 

• Along the boundary, the flow is essentially directed towards that point for all 
times. 
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Such flows can be created by imposing symmetry and other conditions on the initial data. 
For incompressible flows the stagnant point is a hyperbolic point of the velocity field, 
hence the name. 



FIGURE 1. Illustration of hyperbolic flow scenario in two dimensions. 

The scenario is a natural candidate for creating flows with strong gradient growth or 
finite-time blowup, since the fluid is compressed along the boundary. Due to non-linear 
and non-local interactions however, the flow remains hard to control, so a rigourous proof 
of blowup for the 3D Euler equations using hyperbolic flow remains a challenge. The 
crucial issue is to stabilize the scenario up to the singular time. 

One way to make progress in understanding and to gain insight into the hyperbolic 
blowup scenario is to study it in the context of one-dimensional model equations. This was 
begun in Earn, where one-dimensional models for the 2D-Boussinesq and 3D axisymmetric 
Euler equations were introduced and blowup was proven. 

One-dimensional models capturing other aspects of fluid dynamical equations have a 
long-standing tradition, one of the earliest being the celebrated Constantin-Lax-Majda 
model [3]. We refer to the introduction of HU for a more thorough review of known 
one-dimensional model equations, and to [lj for discussion of the aspects relating to the 
hyperbolic flow scenario. 

In this paper, we will study ID models of (jT]) on R with the following two choices of u: 

(3) Ux Hlu , 

(4) u = (-A) - fu; = -c a f \y-x\~^~ a) u}(y,t) dy. 

Jr 

The choice (J3]) leads to a ID analogue of the 2D Euler equation. This model is derived 
simply by restricting the dynamics to the boundary. In section [2] we give a brief heuristic 
argument which works by assuming that u ; is concentrated in a small boundary layer. 

We note that the model (j3J) was mentioned in |T], where it was stated that (J3J) has 
properties analogous to the 2D Euler equation, without giving details. In particular, in 
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[1J a ID model of the 2D Boussinesq equations (an extended version of (J3])) was introduced 
and studied. One of our goals here is to validate the ID model introduced in [I] in a setting 
where comparison with 2D results are available. The fact shown below, that the solutions 
to the model problem (J3]) behave similarly to the full 2D Euler case, provides support to 
the usefulness of the extended version of this model in [I] for getting insight into behavior 
of solutions to 2D Boussinesq system and 3D Euler equation. 

The model defined by (J4j) is called a-patch model and appears in [6], where also a 
viscosity term is present. From the regularity standpoint, the a-patch model is between 
ID Euler u x = Hoj and the Cordoba-Cordoba-Fontelos model u = Hoj (see bed, which 
is an analogue of the SQG equation. These two models differ however from a geometric 
perspective, since the symmetry properties of the Biot-Savart laws are different. For the 
CCF model, the velocity field is odd for even oj, whereas ([4j) is odd for odd oj. It is 
important to choose data with the right symmetry to make u odd, and thus to create a 
stagnant point of the flow at the origin for all times. 

We note that local existence and blowup results for (J4]) were given in [6], where also 
dissipation is allowed. There the authors rely on a suitable Lyapunov function to show 
blowup, whereas we emphasize the more geometric aspects in this paper. That is, we 
will be studying the analogue of the hyperbolic flow scenario for the above ID models 
and show that this leads to natural and intuitive constructions of solutions with strong 
gradient growth and finite-time blowup. 

Another blowup result related to hyperbolic flow was recently proven by A. Kiselev, L. 
Ryzliik, Y. Yao and A. Zlatos [9j and concerns a a-patch model in 2D for small a > 0. 

2. Euler ID model 

2.1. Heuristic derivation. Recall the 2D Euler equations in vorticity form 

OJ t + u ■ Vcd = 0 

where u = V ± (—A) _1 o;. 

We first indicate a simple heuristic motivation for the choice ((2]) (see also [I])- Consider 
the 2D Euler equation in a half-space {x 2 > 0} and denote x = (aq, — x 2 ). The x\- 
component of the velocity (up to a normalization constant) for compactly supported 
vorticity oj is given by 

(5) u t (x,t) = - [ — ^-oj{y,t) dy 

J R 2 \y-x\ 2 

where oj has been extended to {x 2 < 0} by odd reflection (oj(x,t) = —co(x,t)). 

Suppose now that oj is concentrated in a boundary layer of width a > 0 and that 
oj(xi,x 2j t) = oj(x\,t) in this boundary layer. Then a calculation gives 

(6) Mi(aq,0,t) = -2 f log (—— ^ u(yi,t) dy x . 

Jr \ \y i ~ ^lr / 

If we now retain only the singular part of the kernel log ~ — 2 log \z\ and identify 

u with «i, we get (dropping the constants) 
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So a reasonable ID model is 

(7) oo t + uoo x = 0, u x = Pdoo, oo(x,0) = oo 0 (x) (x, t) G M x [0, oo). 
where H is the Hilbert transform, using the convention 

Hoo(x,t) = P.V.f dy. 

J x-y 

For this model, we have the following local well-posedness property: 

Proposition 2.1. Given initial data ooq G 77™((0,1)) with m > 2, there exists T = 
T(||n;o||ir™) > 0 such that the system has a unique classical solution uo G C([0, T\] H 

The proof is standard so we skip it here. 

An alternative argument to motivate Q is to observe that the gradient of the 2D Euler 
velocity is given by a zero-th order operator acting on oo. In one dimension, this leaves 
only the choice u x = cHoo or u x = cut , c being a nonzero constant. So we could also 
consider the model 

(8) oo t + uoo x = 0, u x = —oo, oo(x, 0) = ojq(x) (x, t) G 1 x [0, oo). 

(ED is however not a close analogue of 2D Euler (see Remark 12.51) . 


2.2. Sharp a-priori bounds for gradient growth. We will first prove the global reg¬ 
ularity of the solution to equation ([TJ) by showing that oo x can grow at most with double 
exponential rate in time. Then we will give an example of a smooth solution to ([TJ) where 
such growth of the gradient of oo is achieved, meaning the bound is sharp. 

Due to the Biot-Savart law relating u and oo, the proof of an upper bound for ||u; x (-, t ) ||oo 
is very similar to the proof for the full 2D Euler equations. For the reader’s convenience, 
we give the proof. Recall first the definition of the Holder norm 


\LO\ |c“ = Slip 

\x-y\<i,x^y 


Ms) -u{y) I 

\x — y |“ 


for compactly supported oo. 

We will need an estimate on the Hilbert transform: 


Lemma 2.2. Let 0 < a < 1. Suppose suppioo) C [— D{t),D{t)\ and assume without loss 
of generality that ||o;o||l°° = 1. Then 

IKIloo < C(a) (1 + I log(D(t))| + log(l + IMIc-)) 


Proof. For any 8 > 0, we have 

f My) 


[-D(t),D(t)]\(x—S,x-\-S) ^ V 

Using the oddness of -, we have 


dy 


rD(t) 1 


<C I - dy < C(\ log5| + | log(D(f))|). 
1 s v 


r + u(v) 

P.V. / dy 

Jx—s x-v 


r»r z-\-S 


c—S ^ V 


(u(y) -w(x)) dy 


< C(a)\\uj(x,t)\\c«8 c 


Choosing 8 = min <{ 1, f, we get the desired estimate of 


w. 


□ 


The following Lemma gives an estimate on D(t). 
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Lemma 2.3. Suppose the support of ojq is in [—1,1] and ||u/o||l°° = 1. Then the support 
of co(x, t ) will be inside [•—C exp (Ce ct ), C exp(C'e G4 )], for some universal constant C > 0. 

Proof. Suppose suppo/ = [—D(t),D(t)]. Then for any point x inside of this interval, we 
have 

rD(t) p2D(t) 

|u(x)| < / | log \x — y\\ dy < C / | log |s|| ds < CD(t)(\ log(.D(i))| + 1). 

J-D(t) Jo 

By following the trajectory of the particle at D(t), 

D'(t)<cDmiogm))\ + i)- 

A simple argument using differential inequalities shows that D(t) is always less than z(t), 
where z(t) is the solution of 

z'(t) = Cz(t)(log z(t) + 1), z(0) = min{.D(0), 2}. 

This yields the double-exponential upper bound on D(t). □ 

The following Theorem gives the double exponential upper bound for ui x . 


Theorem 2.4. There is universal constant C such that if oj o is smooth, compactly sup¬ 
ported with suppojQ C [—1,1] and ||cd= 1, 

(9) log(l + ||o/ x ||l°°) < CTog(l + || (^o) x ^L°°)^ Ct (t > 0). 


Proof. We follow the proof in [TO]. Let us denote the flow map corresponding to the 
evolution by Then 

d 


dt 


$t(x)=u($ t (x),t), $ o (a0 = x, 


and 


<9 t |$t(o;) - ®t(y )I 


A 11 11 L°° • 


|$t(®) - $ t (y)| 

After integration, and by Lemma 12.21 and Lemma 12.31 this gives 

for' < l *‘ ( j ) ~*‘ (y)l < /(f), 


\x - y | 


where 


f(t) = exp I C J (1 + exp(Cs) + log(l + Ho/^Hloo)) ds 
This bound also holds for ‘hjT 1 - On the other hand, 


\\u x \\l°° = sup ■ 

x^y 

Which means we have 


\x - V | 


< ||(o/ 0 ) x ||sup 

x^y 


0*0 (y) I 

\x - y | 


(1 + IMloo) < (1 + ||(o/ 0 )x||l°°) exp (C 1 + exp(Cs) + log(l + ||u; x ||z,.x>) ds , 


or 


log(l + 11o/xll^oo) < log(l + ||( 0 / 0)2 1|l°°) + C exp(Ct) + C / (1 + log(l + ||o/ x ||lcx))) ds. 
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So y(t) := log(l + ||cu x ||l°°) satisfies the integral inequality 

y\t) < y( 0) + Ce ct + [ (1 + y(s)) ds 

Jo 

and by the integral form Gronwall’s inequality and some elementary manipulations, we 
arrive at the bound y(t) < C\y(f))e C2t . This yields the desired bound on 1111oo- □ 


Remark 2.5. If we choose our Biot-Savart law to be u x = — u, then from a modification of 
the above proof we get an exponential upper bound for Hco^; - This is different from the 
2D Euler equation, which suggests that (0) is a better analogue of the 2D Euler equation 
than (JED- Moreover the equation (J8l) also has different symmetry properties. 


Next we construct initial data oj$ such that 11(•, t)\\ Laa grows with double-exponential 
rate, proving the sharpness of the a-priori bound (J9|). The hyperbolic flow scenario is 
created in the following way. First, we require that the initial data ojo is odd with respect 
to the origin, and has compact support. By Proposition 12.11 the oddness is easily seen 
to be preserved by the evolution. Consequently, the velocity held (which is also an odd 
function) can be written as 


( 10 ) 

where 

(ii) 


u(x,t) = -X f°° K (-) dy (x > 0), 

«/ o \ y j u 


K(s ) := - log 
s 


s + 1 
s — 1 


Note that the origin is a stagnant point of the how for all times. By taking uq to be 
positive on the right, the direction of the how is towards the origin. More precisely, ujq is 
defined as follows (see Figure [2]): 

• Let cuo be supported on [—1,1], smooth and odd. Choose numbers 0 < xi(0) < 
2x2(0) < 1 such that MxiifS) < 2 ) 2 ( 0 ), where M will be determined later. Require 
that ujq is increasing on [ 0 , 2 q( 0 )], decreasing on [x 2 (0),l] and identically 1 on 
[xi(0),x 2 (0)]. 

Using the earlier notation <L f for the how map associated to (J7j), let 


2)1 (t) := <f>i(xi(0)) 
2)2 (t) := $i(x 2 (0)) 


It is easy to see that the general structure of co 0 will be preserved by the how: For 
hxed t , u>(x,t) will be increasing on [0, Xi(i)], decreasing on [ 2 ) 2 (t), 1] and identically 1 on 
[xi(t), x^ft)]. In fact, since u(x,t ) < 0 for x > 0, xft) and 2)2 (£) will be moving towards 


the origin in time. We will show that the quantity 


3)2 (t) 

2)1 (t) 


increases double exponentially 


in time. This is sufficient to conclude the desired growth of ||u; 2 ,(-, t) H^o 


Theorem 2.6. Assume our initial data is defined as above, then 

lo S exp(CR) (t > 0), 


2)1 (t) 


2)1 ( 0 ) 
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for some positive constant C. As a consequence, 

log ||coL;(•, t) ||L 00 > Ciexp(C 2 1) (■t > 0 ) 


for some C\, C 2 > 0 . 

Theorem 2.4 quickly follows from the following Lemma: 

Lemma 2.7. Suppose 1 > x 2 > 8 x 1 . There are universal constants Co and C 1 so that 


Proof. First observe 


d f x 2 \ 
dt yxi / 




X2 


XL 




d f x 2 \ 

dt yxi J 


x' 2 Xi — x\x 2 u(x2)x 1 — u(x i)X2 X2 



X\ 


u{y) 


dy. 


( u(x 2 ) 
V x 2 


u(x i) \ 
X\ J 


We decompose the integral into 4 pieces which we will estimate separately: 



=/ + // + III + IV. 
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For /, we use 0 < co(y) < 1 and 2x\ < x 2 < 1: 


?2x\ i 

0</< / -log- , 
Jo xi fi - y I 


(xi + y) 


dy + 


= — 3xi log 3 + — 
Xi x 2 

< 3 log 3 + 2 log 2. 


p2xi 

lo 

1 + ^ 


log 


(%2 + y ) 


dy 


2x i lo § , _ Z + ^2 

X 2 


\x2-y\ 

log (l - *1 


V 


X 2 


+ x 2 log ( 1 + — 


X 2 


Using the fact that K(s ) is increasing in [0,1) and decreasing in (1, oo] and that co(y) = 1 
for y G (2xi, |x 2 ) we get 


tX2 


II = 


<2x\ 


K [ — 


y 


^ dy > r 


y 


I2xi 


(2 — x lo g( 3 ))~ 

z y 


(2 — x l°g(3)) log ( 


^2 


Using the positivity of A', 

/•2a;2 ‘ 

///= / A' 


2R 

1/ 


\Ti 


— A' 


-a 


u{y) 


r2x2 


dy>- 


K 


2 X2 


—) w(j/)- dy 

y J y 


> 


4 log 


1 s 


s A 1| 
s - II 


ds > -C. 


We estimate JU in the following way, using that 00 (y) < 1 and < y < 1 for 2^2 < y < 1: 


|/U| = 


'2X2 L 
/■1 


X I — 


< 


y 

y + x 2 


K I — 




y 


dy 


< 


'2X2 L 


A' I — 


K | — 

y 


- dy 

y 


I —log 

I2x 2 V %2 


dy- f -log^I 

.72x2 x i v ~ x i 


= u - u* 


We can compute (i) directly and get 

(i) = — log I + X2 A log(l A x 2 )(l x 2 ) — 2 log(x 2 ) — 3 log(3). 

x 2 1 - x 2 

Similarly, for (ii), we have 

(ii) = — log 1 + Xl + log(xi A 1)(1 — Xi) - 2— log ^' l2 + — - log(2x 2 A Xi)(2x 2 - Xi). 

Note that in the expressions for (i) and (ii), all terms can be bounded by universal 
constants except for —21og(x 2 ) and log(2x2 A xi)(2x 2 — xi). ffowever, using x\ < x 2 , we 
get 

\IV\ < C - 2 log(x 2 ) A log(2x 2 A Xi)(2x 2 - Xi) = C A log ^4 - 

□ 
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The proof of Theorem 12.61 is now completed as follows: choose M > 8 so large such 
that | log (M) — C 0 > 0. We have thus | log ^ j — C 0 > 0. From Lemma 12.71 it follows 

that is growing in time and that we have 


d /x 2 

dt V 


C i x 2 , 

> --log 

2 x\ 



1 


or ^ log (j^J > ~f log (j^j f° r times. This clearly implies that grows double- 
exponentially. 


Remark 2.8. In [10], the Biot-Savart law is decomposed into a main contribution and 
an error term. In our case DO, the main contribution would be 


( 12 ) 



dy. 


If we replace dm hm, then double-exponential growth of can be proven by a straight¬ 
forward argument. In this case, the computation for the estimate in Lemma \ 2. 7| becomes 
much easier. 


3. a-PATCH ID MODEL 

In this section, we consider the ID model equation 

(13) 0J t + ulo x = 0 
with a different Biot-Savart law 

(14) u(x,t) = (-A)- a/2 u(x,t) = -c a f \y - x\~ {1 ~ a) u(y,t) dy, a G ( 0 , 1 ) 

Jr 

For convenience, we will assume the constant c a associated with the fractional Laplacian 
is 1 , and we write 7 = 1 — a. 

This problem has been studied in [ 6 ], where local existence and uniqueness results for 
smooth initial data are proven. From these, we can show that this equation preserves 
oddness and u(0,t) = 0 holds with odd initial datum. For odd data, we can write 

poo 

(15) u(x,t) = - k(x,y)u(y,t ) dy 

Jo 

where k(x, y) = \y — x\ _1 — \y + x| -7 . Note that k(x, y) > 0 for x 7 ^ y G (0, 00 ). 
Following similar ideas as for ID Euler, we specihy our initial data ujq as follows: 

• Pick 0 < 07 ( 0 ), 2 : 2 ( 0 ) with Mx i(0) < £ 2 ( 0 ). Let ojq be smooth, odd, a;o(^) > 0 
for x > 0 and have its support in [—2x2(0), 2 x 2 ( 0 )]. M > 1 is to be chosen below. 
Moreover, let coq be bounded by 1, smoothly increasing in the interval [0,Xi(0)] 
and co 0 = 1 between Xi(0) and x 2 (0). 

As long as the solution remains smooth, the general structure of the solution does not 
change. Let xi(t),x 2 (t) be again the position of the particles starting at xi(0),x 2 (0). 

Theorem 3.1. There exist a choice of xi(0), X 2 ( 0 ), M and a time T > 0 such that the 
smooth solution of (1T31) for the above initial data cannot be continued beyond T. Provided 
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the solution remains smooth on the time interval [0, T), the particle starting at xi(0) 
reaches the origin at time t — T, i.e. 

(16) limxi(t) = 0 . 

t-+T 

In this sense, the solution forms a “shock”. 

Remark 3.2. In j6j, the existence of blowup solutions to (TT3T) is shown using energy 
methods. The advantage is that they are able to include a dissipation term. However, it is 
difficult to see the geometric blowup mechanism clearly using energy methods. Our proof 
for the inviscid case uses the dynamics of the solution and gives a more intuitive picture 
of the blowup, and is easily generalized to other even kernels having the same singular 
behavior. 


In the rest of this section, we will prove Theorem 13.11 So assume that for arbitrary 
choice of xi(0), £ 2 ( 0 ), M, we have a smooth solution u defined for all times t G [0, 00 ). 

First of all, we track the movement of the particle starting at £i(0), which is the 
following Lemma. 

Lemma 3.3. There exists a universal constant M > 2 so that if Mx\(t) < X 2 (t), the 
velocity at x\{t) will satisfy 

(17) u(^i(f),f) < -Cxft) 1 ^, 

for some universal constant C. 


Proof. Let U\ = u(xi(t),t). Since k,oa > 0 on (0, 00 ) 


fX 2 


—Ml > 


k( x uv) dy 


' 2 xi 


= C 7 [~{x 2 + Ti) 1 7 + (x 2 ~ Xi) 1 7 + (3T1) 1 7 - x\ 7 ] 
= c 7 [(3 1-7 - l)xi -7 + (x 2 - Xi) 1-7 - (x 2 + or) 1-7 ] 


= CyX j 


1-7 


(3 1 7 - 1) + ((x 2 - xf 1 7 - (x 2 + xf 1 7 ) 

X, ' 


for some constant c 7 > 0. Note that (3 1 7 — 1) > 0. We can write 


1-7 


T 3 y( x 2 - X1) 1 7 - (x 2 + T1) 1 7 = 

X 1 X-\ 


l- x -l 

X 2 


1-7 


1 + — 

X 2 


1-7 


X 


1-7 


= : '~[zzf (x!/x 2 ). 


x 1 


There exists a constant C > 0 with \ffxi/x 2 )\ < C\xi/x 2 \ for \x\/x 2 \ < 1 / 2 , and so 

-mi > c 7 x }" 7 [(3 1 - 7 - 1) - CM -7 ] 

if Mx\(t) < x 2 (t). Now choose M large enough so that CM _7 is smaller than the number 
|(3 1-7 — 1). “ □ 
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This estimate of velocity field will lead to a blowup in finite time, provided we can show 
Mxi(t) < X2(t). More precisely, 

= u(x i) < 

implying 

®i < c( Xl (oy - ct)$. 

This shows that no later than To := C' _ 1 a;i(0) 7 , the particle x\(t) will reach the origin, 
and the solution cannot be continued smoothly. Note that T 0 does not depend on x 2 (0). 
ft remains therefore to control the motion of x 2 (t), concluding the proof. 

Lemma 3.4. For x 2 ( 0) large enough, Mx\(t ) < x 2 (t) fort e [0, T 0 ). 


Proof. We write u(x 2 (t),t) = u 2 . Observe that the support of u>(-,t) is always contained 
in [—2 x 2 (0), 2 x 2 (0)] because of u(x,t ) < 0 for x > 0. 

Next we find an upper bound on u 2 : 


(18) 

Hence, 

(19) 


\u2(t)\ < 


f 2x2(0) 


'-2X2(0) 


|y — x\ 7 < Cx 2 ( 0) 1 7 . 


To 

x 2 {t)>x 2 (0)- f |u 2 (s)| ds > x 2 (0)(l - (Tr 2 (0)' 7 T 0 ). 
Jo 


Now choose x 2 (0) so large that Mx i(0) < a; 2 (0)(l — Cx 2 ( 0) 7 T 0 ). But then 


Mxi(t) < Mx i(0) < x 2 (0)(l - Ca; 2 (0) 7 T 0 ) < x 2 (t), 


giving the statement of the Lemma. 


□ 
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